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Abstract: We present high precision measurements of polarization 

f-») rotations in the frequency range from 0.1 to 2.5 THz using a polarization 

{S) modulation technique. A motorized stage rotates a polarizer at ~ 80 Hz, 

5_i and the resulting modulation of the polarization is measured by a lock-in 

technique. We achieve an accuracy of 0.05° (900 /zrad) and a precision 
J> of 0.02° (350 /irad) for small rotation angles. A detailed mathematical 

description of the technique is presented, showing its ability to fully 
L^J characterize elliptical polarizations from 0.1 to 2.5 THz. 
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Picosecond (10~ 12 s) timescales are one of the most ubiquitous in condensed matter sys- 
tems. Scattering times of electrons in metals, resonant periods of electrons in semiconduc- 
tors and their nanostructures, vibrational frequencies of molecular crystals, superconducting 
"Cooper pair" decoherence times, lifetimes of collective vibrations in biologically important 
proteins, and even electron transit times in Intel's THz transistor are all phenomena occurring 
in the picosecond range. Such ubiquity makes measurement tools employing terahertz (10 12 
cycles/sec) electromagnetic radiation potentially quite useful. Unfortunately, measurements in 
the terahertz (THz) spectral range have traditionally been challenging to implement as they 
lie in the so-called "terahertz gap" - the range of frequencies above the capabilities of tra- 
ditional electronics, but below that of optical generators and detectors (photonics). In recent 
years, however, a number of dramatic advances in the form of time-domain THz spectroscopy 
(TDTS) have helped span this gap, creating an emerging area of optical and materials research, 
with a broad range of applications (TJE). 

Photoconductive switch based TDTS works by sequential excitation of source and detector 
semiconductor (GaAs in the present work) structures by a femtosecond pulsed laser. Laser 
excitation of a biased Auston switch THz emitter creates an approximately picosecond electro- 
magnetic pulse (with THz frequency Fourier components) that propagates through space and 



interacts with a sample under test. Detection is accomplished via another photoexcited Aus- 
ton switch, which is biased across its electrodes by the transient THz electric field. With these 
developments, THz spectroscopy has become a tremendous growth field [3|, finding potential 
use in a multitude of areas including materials characterization for solid-state devices J4] 13, 
optimization of the electromagnetic response of novel coatings [6|, probes of superconductor 
properties QUO, security applications for explosives and biohazard detection |9], the detec- 
tion of protein conformational changes [10], and non-invasive structural and medical imaging 

Despite the advances in THz spectroscopy, a number of challenges remain. For example, 
highly accurate measurement of polarization states in the THz range has proven to be difficult 
lfl4l . Historically the highest achievable precisions in the THz or millimeter wave range have 
been approximately 1°, measured using rotatable wire grid polarizers, with higher precisions 
only being achieved with extremely specialized setups 1 15 1. This is nowhere near the sub-/irad 
resolution that is possible in the near infrared and visible range |[T6l[T7l . A number of groups 
have developed polarization sensitive switches [18, 19 1 fabricated in a multi-pole configuration, 
where a single device is sensitive to two orthogonal electric field polarizations. This simultane- 
ous detection of both electric field polarizations is advantageous, as the phase sensitive nature 
of TDTS means that unlike conventional polarimetry, only two orthogonal directions have to 
be measured to resolve the complete Jones matrix for a polarized wave. 

Although multi-pole devices lfl8l [191 and calibrated wire-grid polarizer measurements GUI 
are capable of angular precision of 0.2-0.3°, one might hope for even higher angular precision, 
as many material systems exhibit very small Kerr or Faraday rotations that are of fundamental 
interest for condensed matter physics |[T7ll2"Tll2"2"ll2"3~ll24l . The possibility of even higher preci- 
sion THz polarimetry would also open the door for entirely new spectroscopic techniques such 
as THz ellipsometry 1251 . Currently the best THz ellipsometers (based on continuous wave 
sources) cannot take reliable data below ~ 3 THz ||26l . 

With these ends in mind we have investigated and applied a polarization modulation tech- 
nique using a fast rotator in combination with photoconductive switch based time-domain THz 
spectroscopy . This technique, originally developed for use with continuous wave gas lasers and 
Fourier transform infrared spectroscopy [27 28 1, shows a distinct advantages when applied to 
TDTS, namely that the phase sensitivity of TDTS allows the unambiguous measurement of the 
Jones vector polarization of a THz wave. The technique has recently been applied to TDTS 
using electro-optic detection on topological insulators [29, 30 1. We find that when used with 
photoconductive switch TDTS it is capable of an unprecedentedly high angular precision of 
better than 0.02°, with an accuracy of 0.05° in the frequency range from 0.1 to 1.25 THz, 
a sensitivity that makes entirely new classes of measurements possible. In this paper we lay 
out the mathematical analysis of a polarization modulation measurement and describe initial 
experimental characterizations of this technique. 

1. Theory 

In this technique, a rotating polarizer is used to modulate the polarization state of the terahertz 
light. In a manner to be described below, it allows for direct detection of the vertical (x) and 
horizontal (y) components of the electric field after passing through a sample using the in-phase 
and out-of-phase channels of a lock-in amplifier, respectively. The advantage of this technique is 
three-fold. First, measurement time is reduced as a single measurement determines both electric 
field polarization states, compared to the two measurements required with static polarizers. 
Second, because the two polarization channels are measured simultaneously, time dependent 
common mode noise may be ratioed out, improving the signal of measurements such as Kerr 
and Faraday spectroscopies. Third, as the experiment is only sensitive to the modulation of the 



polarization state, many detrimental effects due to the finite extinction ratio of the analyzing 
polarizers cancel out, improving the accuracy of the measurement. 

We first describe the mathematical analysis in detail. As usual, a time -domain pulse E x (t) 
can be written as a superposition of Fourier components 

E x (t) = -^L f d(0e im E x {(0) (1) 

\/2k J-<*> 




Here, E x (t) is the purely real electric field propagating through the system, and E x (to) de- 
scribes the amplitude of the complex Fourier components that comprise E x (t). The mathemati- 
cal description of the pulse as it propagate through the system is easiest to perform in terms of 
these Fourier components. The propagation of the Fourier components as they travel through 
the system can be described with the Jones matrix formalism. In this formalism a matrix with 
(in general) complex elements acts on the Fourier component with electric field amplitude E° 
to produce a final polarization state E f : 

M xx (a) M xy {(0) \ ( E°(ffl) \ 

M yX {(0) Myy{(0) )\E*{<0) ) {Z) 

The advantage of this formalism is that a complicated series of optical elements can be rep- 
resented as a simple set of matrix multiplications M on the original electric field Fourier com- 
ponent. 

The Jones matrix for a polarizer with its transmission axis oriented at an angle © with respect 
to the x axis is given by: 

/ cos 2 (0) cos(0)sin(0)\ 

^cos(0)sin(0) sin 2 (0) J K ' 

Here, a polarizer with = polarizes light in the x direction and is denoted P x , and a polar- 
izer with © = n/2 is called P y . For a rotating polarizer with angular velocity SI, the angle is 
simply replaced by the factor Sit to obtain the polarizer matrix as a function of time. 

In general, the goal of ellipsometry is to determine the values of the complex elements of the 
Jones transfer matrix for a sample. The changes in amplitude and phase these elements produce 
in the electric field can then be correlated with physical properties of the sample, such as the 
conductivity. The complex transfer matrix T for a sample is represented as: 

V ' \ty X (0}) tyy(CO)) 

In the setup described here, a rotating polarizer is used to modulate the polarization of a tera- 
hertz waveform after it has passed through the sample (Fig. [Til. Analysis of the time dependence 
of the waveform can be used to extract the frequency dependent components of the T matrix. 

Using the Jones matrices to describe the frequency domain Fourier components of the wave- 
form as they travel through the sample and rotator, the result is: 



E® [txx cos 2 (Sit) + t yx cos (Sit) sin (Sit)] 
E® [t xx cos (Sit) sin (Sit) + t yx sin 2 (Sit)] 

+£° [t xy cos 2 (Sit) + t yy cos (Sit) sin (Sit)] \ 

+E° [txy cos (Sit) sin (Sit) + t yy sin 2 (Sit)] ) '""' 



These electric fields have a mixed dependence on the components of the T matrix. The 
analysis can be greatly simplified by adding vertical polarizers before and after the rotator and 
sample, respectively, (Fig.fTT) to isolate t xx and t yx : 



P x Pn,TP x - 



?0 



f El (t xx cos 2 (Qt) +t yx cos (Qt) sin (Qf))\ 
) 



E' 



| {t xx ( 1 + cos (2Qf )) + f yx sin (2Qf )) j (6) 



or alternatively adding horizontal polarizers to isolate t xy and f v> ,: 



PyPbt-T-P, , , 



v ^ (?„ sin (2flf ) + f„. ( 1 - cos (2Qf ))) 



(7) 



Note that photoconductive antennas themselves have polarization sensitivity and naturally 
the best results will be obtained if they are oriented along the direction defined by the static 
polarizers. 

To understand the detection of the waveforms by the photoconductive antennas, there are 
two relevant times to keep track of in the experiment. First, there is the time associated with the 
position of the delay stage (Fig. [TJ, determining the small time window of the terahertz pulse 
that is being measured. We denote this time as t„. This is connected to the position of the delay 
stage L by t p — 2L/c. Second, there is the time associated with the rotation of the polarizer, 
denoted as t r . This is the time that is referred to in Eqs. (J6l) and (FT}. We can now rewrite Eq. (fT) 
for the purely real electric field measured at the detector using this notation: 

E{(tp,t r ) = 4= r dcoe im "Ef(co,t r ) (8) 

V 2.K V-oo 

where E x (co,t r ) is defined in Eq. (roll. In the lock-in amplifier this signal is then mixed with 
the second harmonic of Q. and integrated for an amount of time determined by the lock-in time 
constant, T. The signal of the in-phase lock-in channel is given by 

Sx{t p ) = - f dt r cos(2&t)E{(t p ,t r ) 

T JO 

= —= [ T dt r cos(2Q.t) [ dco'e im ''"E{(co',t r ) (9) 

TV2ft ^o J—*> 

Taking the Fourier transform of Sx (t p ) allows the Fourier components to be extracted: 

Sx(C0) = — f dt r cos{2£lt) I dt p e~ iat P [ d(o'e im ' t i'E{((o',t r ) 
2kt Jo J—*> J-°° 

1 f T 
= - / dt r cos{2D.t)E{(co,t r ) (10) 

i Jo 



Substituting in for Ex(co,t r ) from Eq. (|6ji: 



+ £ " (fi)) 2 f -" (fi)) (sin (20f r ))) cos (2Qr r ) 



T V 4 Jo 



dt r (2 cos (2i2f r ) +cos (4£2f r ) + 1) 



■^WW / T ^sin(4n f / 
4 Jo 

= ~tf (a>)E° ((»)?«(<») (11) 

where /?o (fi>) parameterizes the combined responsivity of the antenna and lock-in amplifier. In 
the y channel of the lock-in the signal is mixed with sin (2£it r ), and similarly gives: 

S Y (co) = ^R (co)E x (co)t yx (co) (12) 

How do these measured values relate to the quantities we wish to measure, namely the elec- 
tric fields after the sample? An analysis of the effect of the sample on vertically polarized light 
without the rotator shows that the electric fields after the sample are: 




Comparing the results of Eqs. (Hi and ( 12 1 with the fields found in Eq. ( 13 I, it is apparent 



that the rotator measurement is equivalent to direct measurement of the x and y components 



of the electric field after the sample simultaneously in the two lock-in channels. Eqs. (Ill and 



( 12 1 represent the principle result of this analysis: that a single measurement with the rotator 
technique is exactly equivalent to two measurements with a standard static polarizer setup. 

In a real system, the non-ideality of the components must be taken into account. For exam- 
ple, perfectly aligned off-axis parabolic mirrors produce no rotation of the polarization state, 
but in a real system perfect alignment is impossible to achieve. In this regard placement of op- 
tical elements in the setup is crucial. In Fig. [T] off- axis parabolic mirrors (OAPs) 1 and 2 can 
introduce a rotation 0i2 of the initial polarization state of the electric field from the emitter. To 
eliminate measurement errors introduced by this rotation, it is critical to place the first static 
polarizer after OAPs 1 and 2, which essentially redefines the initial vertical electric field for the 
measurement at this point. 

For OAPs 3 and 4 the same problem exists. Consider the situation explicitly if polarizer P2 is 
placed not immediately before OAP 3, but instead after OAP 4. The total rotation 034 produced 
from small misalignments in OAPs 3 and 4 can be represented by the rotation matrix: 

COS (0 34 ) Sin (034) \ (14) 

-Sin(0 34 ) COS ( 034 ) ) 

The x polarized light detected in the two channels of the lock-in becomes a mixture of the 
components of the T matrix: 



R E° 

S X = — — {fxx COS ( 034 ) - tyx Sin (034 ) ) 

R E° 

Sy = -^ L (fvxcos(03 4 )+f rx sin(03 4 )) (15) 

As before, placing a polarizer between the rotator and OAP 3 simplifies the result: 

S x = —^t xx cos (03 4 ) , S Y = — —t yx cos (0 34 ) (16) 

Now each channel acquires the same rotation factor of cos ($34), which merely acts to scale 
the overall magnitude of the electric field. There is no mixing of the polarization states in the 
lock-in channels. 

The finite extinction ratio of the polarizers in the system is also a factor in the overall system 
performance. The Jones matrix for an imperfect polarizer is [31 J: 






cos 2 (0) + 77 sin 2 (0) (l - tj)cos (0) sin (0) 
v (l - 77) cos (0) sin (0) 77 cos 2 (0) + sin 2 (0) 

where tj is the extinction ratio of the polarizer, defined as the ratio of the transmitted electric 
field when the polarizer transmission axis is perpendicular and parallel to the electric field direc- 
tion, 77 = E'"" 15 / El ra " s . Accounting for the finite extinction ratio (rj r ) of the rotating polarizer 
in the Jones matrix analysis, the measured signals become: 

Sx = ^t xx cos (034) (1 - ri r ) , S Y = ^f^cos (0 34 ) (1 - r] r ) (18) 

The finite extinction ratio of the rotating polarizer acts only as a simple scaling of the am- 
plitude of the measured electric field. This highlights a distinct advantage of the modulation 
technique over the standard static polarizer measurements: only modulated signals are meas- 
ured by the lock-in, and non-ideality of the rotating polarizer only produces a small effect on the 
amplitude of the measured electric fields. Commercially available polarizers in this frequency 
range typically have extinction ratios of tj ~ 1/50, resulting in only 2% difference between the 
ideal and real signals. The best figure of merit for the rotating polarizer becomes its stability 
under rotation and its robustness to high rotation rates to reduce 1 // noise, not the extinction 
ratio. 

Finally, the effects of the finite extinction ratios of the static polarizers must be considered. 
This changes the form of the vertical polarization matrix, P x , to 

Finite extinction ratios of the static polarizers mix different T matrix components in the 
measured signals. Since the rotator projects the electric field onto the second polarizer at large 
angles, a non-trivial electric field leakage results and the elements of the T matrix are mixed 
in the two detection channels. The obvious solution is to increase the quality of polarizer, but 
while polarizers with significantly higher extinction ratios in the terahertz range have been 
demonstrated, they are still not at the stage of commercial production l32l . As shown below 
the solution to this problem relies on the final analysis performed on the collected time domain 
waveforms. 

As described above, in order to obtain the electric field components as a function of fre- 
quency, the measured time domain waveform is Fourier transformed to obtain the frequency 



response. For many measurements performed in this system, the quantity of interest is the ro- 
tation angle introduced by the sample, which can be found by simply taking the ratio of the 



Fourier transforms of S y and S x . For the case of vertical PI and P2 polarizers (Eq. (Hi) 

S y {(o) E° x {co)t yx {(o) tyxica) 



(20) 



S x (co) E° x (©)&*((») fxx((0) 
If the polarizers are assumed to be ideal, simple cases such as rotation can be explicitly 



solved. The T matrix for a rotation is given by Eq. ( 14 1, which allows Eq. ( 20 1 to take a 
simple form 

t yx {co) sin [0(a))] 



t xx {(0) cos [0 (to)] 
which easily gives the rotation angle 



arctanl'^gf )=*,«>, (22, 

\COS[<J)((0)]J 



However a finite extinction ratio of the polarizers complicates Eq. pO) . If the finite extinction 
ratio of the static polarizers in the P x geometry is taken into account, it is impossible to solve Eq. 



(22 1 analytically for the rotation. The situation is particularly complicated for the case where 
small misalignments in off-axis parabolic mirrors given finite angles 0i2 and ^34. However, we 
can perform a Taylor expansion in 77 up to second order for a simple rotation, giving 



Sy (tO) 



arctan y = (to) - - 



S x {(0)J rK ' \ E°-r-£°cot(fr 2 ) 

E° \ 

+ F o^ F o t 72 N + tan (^4) ri (23) 

E x '+E\! tan (<j) n ) J 

This expansion shows that a single measurement is in general not sufficient to determine the 
rotation angle. Even with the setup optimized by proper polarizer placement to eliminate first 
order rotation effects due to the off-axis parabolic mirrors ($12 and ^34) (as discussed above), 
the finite extinction ratio of the polarizers can still cause systematic error. Note however, that 



the second term in Eq. (23 1 does not depend on the rotation angle 0. In this regard, two meth- 
ods can be used to eliminate the systematic error. If the sample itself can be used to give a 
reference angle of rotation (say at a particular temperature or magnetic field), this reference 
can be subtracted to eliminate the constant term in Eq. ( f23) and obtain the relative change in 
rotation angle. An example could be comparing a rotation measurement above and below T c in 
a superconductor for evidence of time reversal symmetry breaking [ 17 1 or of a magnetic ma- 
terial above and below a magnetic ordering temperature. Alternatively, if an accurate reference 
sample can be used, the sample and reference angle scans can be subtracted, and systematic 
errors of this kind can be eliminated. Also note that in general the contribution of E x will be 
small, as in the typical experimental geometry the vertically oriented photoconductive switches 
give an electric field which is primarily oriented along the x direction. 

Up to this point, a technique has been developed and shown to be capable of measuring 
simple rotations. What are the ultimate capabilities of this technique, and how do they connect 
to full terahertz ellipsometry measurements? To answer this question, we must first define what 
is measured in a standard ellipsometric measurement. In standard ellipsometry, light with a 
known polarization is focused onto a sample, and the reflected or transmitted light is collected 



by polarization sensitive detectors. The amplitude of the light in the x and y axes is determined 
as a function of the frequency from this measurement. 



Ey{a>) 



■■ tan[y/(a))]e 



jA(ffl) 



(24) 



E x (a) 

where tany/ = \E y /E x \ and A represents the temporal phase difference between the two field 
components [ 33 1 . The broadband polarization modulation measurement presented here gives 
E x (ffl) and E y (o), so both y/ (©) and A (©) can be calculated and a full frequency dependent 
ellipsometry measurement is possible. 

Next we can ask, how do these measured complex field components connect to rotations and 
the ellipticity? In general, fully polarized states of monochromatic transverse electromagnetic 
waves trace out an ellipse in time with major (a) and minor (b) axes, each at an angle 9 with the 
x and y axes in the laboratory frame, respectively 1 33 ]. The ellipticity is defined by the major 
and minor axes as tan(e) = -. For example, for circularly polarized light a — b and e = ±45° 
with + for right circular polarization, and - for left. Linearly polarized light is characterized by 
b = and e = 0. All other polarization states are called elliptical. For the simple polarizer in 
Eq. (|22|, the ellipticity £ = 0, and the angle 9 corresponds exactly to the angle of rotation (p. 

To perform a complete analysis, it is important to note |34| that the ellipse the electric field 
traces in time is related to the measured ellipsometric quantities tan y/ and A by 



tan(20) = 
sin(2e) = 



2 tan y/ cos A 



1 



tan 2 y/ 



2 tan l// sin A 



1 



tan 2 y/ 
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To start the analysis, we take arctanp defined in Eq. (24 1, which can then be expanded 
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First taking the real part of Eq. ( |27| ) we find 
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Inserting Eq. (|25| into Eq. (|2"8j) gives 



Re [arctanp] = Re 
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which is simply the rotation angle (j) that was found for the simple polarizer case. Now we find 



the imaginary part of Eq. ( 27 1 
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Fig. 1. Experimental setup for the fast rotator measurement. A linearly polarized terahertz 
waveform passes through a sample, introducing ellipticity. The rotating polarizer modulates 
the polarization, and the modulation is detected by the lock-in amplifier. 
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Using Eq. ( 26 1 we can simplify: 
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(30) 



Im[arctanp] = - arctanh (sin (2e)) 



(31) 



The Im(arctanp) can be used to determine the the ellipticity e. This analysis shows that 
the complex angle extracted from arctanp contains all the information necessary to determine 
the frequency dependent elliptical polarization state of the final terahertz electric field. The real 
part is a direct measurement of the rotation angle of the light, and the imaginary part determines 
the ellipticity of the final electric field. In fact, this analysis does not rely on the polarization 
modulation technique, but is valid for the static measurement described by Eq. ( p"3j ) as well. 
The advantage of the polarization modulation technique is that the full characterization of the 
frequency dependent electric field ellipse can be extracted from a single measurement. This 
shows the power of measuring the time dependent electric field and collecting the full amplitude 
and phase information of the waveform. 



2. Experiment 

The experimental setup is shown in Fig. [T] The generation and detection of the terahertz wave- 
form is accomplished by a standard time-domain terahertz spectrometer using photoconductive 
antennas. A ~ 60 fs pulse from a KM Labs TkSapph oscillator with an 80 MHz repetition rate 



is split into two paths via a beamsplitter. The first pulse strikes a biased photoconductive an- 
tenna, generating a terahertz waveform with a spectral bandwidth from 0. 1 to 3 THz. The OAPs 
are laid out in an 8/ geometry. OAPs 1 and 2 are used to collimate and focus the waveform. A 
wire grid polarizer (PI) after OAP 2 defines the initial vertical polarization state of the system 
as the terahertz light is focused onto the sample. The vertically polarized light passes through 
the sample, where it becomes elliptically polarized according to the sample transfer (T) ma- 
trix. The elliptically polarized terahertz light passes through a rotating wire-grid polarizer (RP), 
made by QMC Instruments, which modulates the elliptical polarization at a frequency £2. A 
second vertical polarizer (P2) projects the rotated light back to the vertical (x) axis, and a pair 
of off-axis parabolic mirrors recollimate and focus the terahertz light onto the detector antenna. 
The second pulse created by the beamsplitter is used for detection of the waveform. An addi- 
tional vertical polarizer (P3) is used to account for any rotations produced by OAPs 3 and 4, as 
the antenna has a small, but finite response for the horizontal polarization, which would result 
in mixing of the elements of the T matrix in the detection as shown above. The terahertz elec- 
tric field at the antenna produces a small AC current, and a transimpedance amplifier is used to 
convert this to a voltage. The voltage signal is detected by an SRS 830 lock-in amplifier, where 
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Fig. 2. System resolution characterization, a) Measurement of test polarizer angles from 
20° to 80° . b) Comparison of the accuracy and precision for small and large test polar- 
izer angles, c) Ellipticity angle e for two polarizer angles, d) A test polarizer angle of 6°, 
showing 0.05° accuracy and 0.02° precision up to 1.25 THz. 



it is mixed with the second harmonic of the rotation frequency £2. The rotation frequencies used 
here (2D. between 50 Hz and 80 Hz) are chosen to be high enough to significantly reduce 1 // 
noise, but small enough to minimize the possibility of damaging the polarizer. 

The rotator was built by the Instrument Development Group in the Johns Hopkins Depart- 
ment of Physics and Astronomy. The polarizer is held by a 2" aluminum cylinder that can be 
rotated at high speeds using Bearing Works silicon nitride ball bearings. This is connected via a 
belt drive to a high speed AC motor, the Faulhaber Minimotor 4490, which has a variable speed 
controller that can be used to rotate the polarizer up to 2D = 200 Hz (6000 rpm). To generate 
a trigger signal for the lock-in, two small holes are drilled in the rotating cylinder, and an LED 
and photodiode are placed on opposite sides of the cylinder. When a hole passes the LED, it 
triggers the photodiode on the opposite side of the cylinder. This produces pulses at twice the 
rotation frequency, giving the 2D. signal needed for the lock-in detection. All this is held in 
a metal casing that can be securely attached to the optical table to reduce detrimental effects 
caused by motor vibration. 

To characterize system performance, a wire-grid polarizer in a static rotation mount was used 
as a reference sample. To demonstrate the versatility of the technique, a number of polarizer 
angles were measured with the results shown in Fig. [2] For large rotation angles, the precision 
and accuracy of the measurement is reduced, as the comparison of the measurements for rota- 
tion angles of 20° and 80° (Fig.[2b) shows. This seeming loss of precision and accuracy has its 
origins in the non-ideality of the polarizer as a test sample. The polarizer does not simply rotate 
the polarization, but rather projects the polarization with an angle dependent (cos 2 0) attenu- 
ation. For large angles, the sample polarizer significantly reduces the amplitude of the initial 
terahertz electric field, thereby decreasing the achievable signal to noise. The finite extinction 
ratio of the polarizers also becomes more important for larger polarizer angles, as the ampli- 
tude of the orthogonally polarized field transmitting through the polarizer grows. Additionally, 
at these large angles the amplitude of the orthogonal transmitted electric field becomes highly 
frequency dependent. These two effects are shown in detail in Fig.[2h At 20°, the measured 
angle is constant to within 0.1° over the range from 0.1-1.75 THz, with a precision of ~ 0.05°. 
At 80°, the angle is only constant to 2° within the same frequency range, and the precision of 
the measurement is reduced to ~ 0.5°. The ideal test sample would be a broadband half wave- 
plate that could rotate without attenuation. However, broadband terahertz waveplates are at the 
developmental stage [35 1, so a wire grid polarizer serves as the simplest available test sample. 

As demonstrated above, the accuracy and precision of the technique improve significantly for 
smaller rotation angles of the test polarizer. Fig. |2j: shows a measurement for a polarizer angle 
of 6°. For this small angle, within the range of 0.1-1.25 THz the measured angle is constant 
to within 0.05°, and from the signal to noise ratio, we estimate that changes in the angle as 
small as 0.02° (350 jtirad) can be resolved. The main limitation in measuring the accuracy of 
the technique is the difficulty in setting the angle of the test polarizer to the required level of 
accuracy. At this point we can measure THz polarization steps much more precisely that we 
can generate them. The precision of the system is thus far only limited by the averaging time. 
The scans shown in Fig. [2] each account for 20 minutes of measurement, which is reasonable 
for many experiments. 

To test the technique on a real sample, the birefringent response of a piece of X-cut sapphire 
was measured. The sample was placed with the ordinary and extraordinary axes at 45° with 
respect to the initial vertical (x) light polarization. This crystal orientation projects the electric 
field of the terahertz waveform equally on the two crystal axes, producing a phase delay be- 
tween the two components and changing the polarization from linear to elliptical. For a purely 
birefringent material the T matrix takes the form 
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Fig. 3. a) Single measurement of the birefringence of sapphire with the electric field inci- 
dent on the sample at a 45° angle to the two principle axes, b) Calculation of the birefrin- 
gence (n e — n ) for the ordinary and extraordinary axes for a static polarizer measurement 
(requiring 2 measurements plus a reference) and a polarization modulation measurement 
(requiring 1 measurement). 
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where a is the angle between the extraordinary axis and laboratory x axis, and (p x and (p y are 
the additional phases associated with the electric field traveling along the extraordinary and 
ordinary axes, respectively. Here, with a = 45°, Eq. ( 32 1 becomes 






For the birefringent response, introducing A<p = <p T — qiy, we find that 
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meaning the phase difference between the two axes is 
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In Fig. [3^, we show a time-domain trace of the measured electric field amplitudes as a func- 
tion of time, measured directly on the in-phase and out-of-phase channels of the lock-in am- 
plifier. It shows quite readily the conversion of the initial linearly polarized THz pulse into 
an elliptically polarized state. From the Fourier transforms of the data and the phase differ- 
ence between field components, the frequency dependent difference in the index of refraction 
between the ordinary and extraordinary axes can easily be computed using An = cA(j)/2ndf, 
where d is the thickness of the sapphire. In Fig. [3j? we compare the measured birefringence 
of the sapphire taken with the conventional technique using static wire grid polarizers and the 
rotator technique. For the static wire grid polarizer technique, the ordinary axis of the sapphire 
is oriented along the x axis. A wire grid polarizer oriented at 45° is placed before the sample, 
defining the incoming electric field polarization at 45° to the ordinary and extraordinary axes. 
An analyzing polarizer is placed after the sample, and two measurements are done at ±45°. A 
quick mathematical analysis with the tools described above shows that the two measurements 
then give 

S +A5 o = e^+e'fy (36) 

thus taking the ratios of these two measurements yields the same result for A<p as the rotator 
measurement. 

Fig. [3b shows that the two measurements produce slightly different results. The difference 
in the two measurements is ascribed to the imperfect nature of the polarizers. For the fast ro- 
tator experiment, to second order the effects of the finite polarizer extinction ratios drop out. 



In the static experiment, analysis using the imperfect polarizer (Eq. 17 1 changes Eq. (36 1 such 
that e* — > ( 1 + 7] ) 2 e* and e"k- — > ( 1 — j]) 2 e ! ^ . The robustness of the polarization modulation 
technique actually produces a measurement that is less prone to error than that with the static 
polarizers. This simple measurement demonstrates the power of the technique: a single polar- 
ization modulation measurement produces the same information as three measurements with a 
static polarizer setup, and gives much improved accuracy and precision. 

In conclusion, we have demonstrated high precision measurement of polarization states using 
time domain terahertz spectroscopy. We can resolve angular rotations with an accuracy of ~ 
0.05° and a precision of ~ 0.02°. Additionally, a number of practical mathematical results 
useful in the analysis of phase sensitive polarization measurements have been presented. We 
believe that this technique will have a wide applicability to a number of materials systems at 



the forefront of condensed matter physics, such as high-T c superconductors, quantum magnets, 
and topological insulators. 
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